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Abstract — The present paper deals with an inventory 
management for deteriorating items with instantaneous supply 
and alternating demand rates. A two parameters weibull 
distribution is taken to represent the time to deterioration and 
the demand is changing at a known and at a random point of 
time in the fixed production cycle. Finally, two numerical 
examples have been given for such demand fluctuations to 
illustrate the proposed model. 
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I. Introduction 

The economic order quantity model is the oldest 
inventory management model. In 1915, Harris[l] has 
discussed this model. Ostenyoung et al[2], in one of their 
research publications have reported that this model widely 
used in industry, although the assumptions necessary to justify 
its use such as constant demand rate, known inventory holding 
and set up costs, no shortages allowed etc. are rarely met. 

In many inventory models, the effect of deterioration is 
very important. Ghare and Schrader[3] developed a model 
assuming a constant rate of deterioration. The assumption of 
constant deterioration rate was relaxed by Covert and 
Philip [4] who used a two parameter weibull distribution to 
represent the distribution of the time to deterioration. The 
model was further generalized by Philip [5] by taking a three 
parameter weibull distribution. Misra[6] also adopted a two 
parameter weibull distribution deterioration to develop an 
inventory model eith finite rate of replenishment. In this 
connection, the works done by Jalan et al[7], Dixit and 
Shah[8], Gupta and Agarwal[9], Cohen[10], Hwang[ll] are 
note worthy. 

An extensive research work has already been done by many 
researchers on inventory model by assuming a constant 
demand rate. Moreover the time dependent demand rates have 
been investigated by various researchers like Donaldson[12], 
Silver[13], Ritchie[14], Datta and Pal[ 15], Mandal and 
Pal[ 16] etc. But Dhandra and Prasad[17] have discussed an 
inventory model assuming demand rates for the inventory 
changing at a known and at a random point of time in the fixed 
production cycle. Pal and Mandal[18] have extended this 
work to a situation where the inventory is deteriorating at a 
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constant rate. Such fluctuation in the demand rate is quite 
realistic and it may occur due to any of the following reasons 

(i) . a particular item becoming obsolete day by day and 
non- 
availability of its spare parts 

(ii) . Due to higher or lower purchasing cost of items. 

In the present paper, we derive an inventory model of 
items that deteriorate at a weibull distributed rate, assuming 
the alternating demand rates. Moreover the model has been 
developed for both the deterministic and stochastic situations. 
Some particular cases have been discussed and some 
numerical examples have been provided for such both 
demand pattern to illustrate the present model. 

II. Preliminaries 

Assumptions: 

(i) Lead time is negligible. 

(ii) Replenishment size is constant and its rate is 
infinite. 

(iii) Shortages are not allowed. 

(iv) A deteriorating item is neither repaired nor 
replaced during a given cycle. 

(v) Length of each production cycle is constant. 

Notations: 

T : Length of one cycle. 

I(t) : Inventory level at time t( > 0). 

R(t) : Demand rate at time t( > 0). 

Y : Time epoch at which the demand rate changes. 

K : Fixed ordering cost of inventory, 
c : Cost of unit item, 
h : Holding cost per unit per unit time. 

0(t) : The deterioration rate function for a two parameter 
weibull distribution. 


III. Mathematical formulation and solutions 

Let Q be the total amount of inventory produced at the 
beginning of each period. Since we have assumed an 
alternating demand rate for the cycle time T, we can consider 
the demand rate R(t) in the following form: 
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R(t >= {t 


0 <t<z 
z<t<T 


where z=min(Y,T) is the epoch at which the demand rate 
changes. 

The deterioration rate 0(t) is taken as 

0(t) = H p t P ~' , 0<a«l, P > l,t > 0. 

Hence the differential equations governing the system in the 
interval (0,T) are the following 


I(t) = Q(1 -at P )- A [t- 


a/? 

JV\ 


t p+\ 


],0<t<z 


(4) 


and 

a 


I(t) = Q(1 - a t P ) - ( A l - A 2 )[ z(l - a t P ) + 


P + 1 


^ +1 ] 


[ t ■ 


a/? 

JV\ 


,p+\ 


], z<t<T 


(5) 


dl(f) 

dt 


+ aP t p 1 1(t) = - A, , 0 < t < z 


0 ) 


Using the condition I(T) = 0 and neglecting the terms of 
0(a“ ), we have 


dl(t) ^ ^ 

and — — +ap t p I(t) = - A 2 ,z<t<T (2) 
dt 


Q = a-^)[z+ J?— Z P +i } + A,[T+ ~^—T P+l 


p+l 


p+l 


] ( 6 ) 


The boundary conditions are 1(0) = Q and I(T) = 0 (3) 

2 

When 0<oi«l, we ignore the terms of 0( OC ), and use the 
condition 1(0) = Q , then the solutions of the equations (1) and 
(2) are respectively the following 

The average number of units AI(T) in the inventory during the period (0,T) is 

AI(T) = — [ T I (t)dt = — [ \~ I (t)dt + \ T I it)dt] 

T Jo y Jo Jz 

Substituting the values of I(t) in the above integrals and then eliminating Q using equation(6) and neglecting higher order terms 

of a we find the following 


1 , , z 2 aft 

AI(T) = — + 

T - 2 (P + \)(P + 2) 


V P + 2 


} + A 2 { ■ 


ap 


-P+2 


C(T, z) 


2 (P + lXP + 2) 

Therefore, the average total cost per unit time in (0,T) is 

K + cQ h r / , 0x2,, 2a/? 


}] (7) 


+ — 1(A 1 -A 2 )z 2 { l + 


(P + 1XP + 2) 


z"j 


A 2 T 2 { 1 


2ap 


(/? + !)(/? + 2) 


T p }] 


— + - [(A l -A 2 )z{ l + z p } + A 2 T{ 1+ 

T T p + l ~ p + l 


T P }] 


^[{\-A 2 )z 2 { l 


2 ap 


(/? + !)(/? + 2) 


z p }+ A 2 T 2 {\ 


2a P 


(P + l)(P + 2) 


T p }] (8) 


We now minimize the above total average cost per unit time under the following two cases 

(i). z is a known point of time 
and (ii). z is a random point of time. 

Case-1: z is a known point of time. 

(a) When z > T. 

In this case, the expression for C(T, z) reduces to the following 

2a/? 


K a hA, 

C(T) = — +cA.[ 1 + T p ]+ — l T[ 1 + - 

T 1 p + l 2 (p + l)(P + 2) 


T p ] 


which is the cost function of the economic order quantity inventory model for deteriorating items. 
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If there be no deterioration i.e. a = 0, the above expression becomes 

K „ h / L 

C(T) =-+c/L + — - T 

T 2 

which is similar to the equation obtained by Harris and Wilson. 

(b) When z<T. 

We consider z = pT, 0< p < 1. 

Thus from expression(7), we find 


afi/Ll 1 


-T p+2 }+ A,{— + - 


a(5 


1 T 2 

AI(T) = — [ ( /lj - /t, ) /j 2 { — + 

T “ 2 (/? + 1)(/? + 2) 2 (/? + l)(y? + 2) 

Therefore the equation (8) reduces to the following 
QT) = J* c - ) 1 + + ^T{ 1 + 1] 

+ ± [( r,-^)^rPi + 2a ^ ^ 

27 1 " (J3+1XJ3 + 2) 


T p+ 2 }] 


T p } + JL T 2 \ 1 + 


lap 


00 + 1)00 + 2) 


r^}] (9) 


This gives 


dC(T) 

For minimum average cost, the necessary condition is = 0 

dT 

ha/3 -T p+2 [\ l u l}+2 +x 2 {\-/ +2 )] + ^_ rnv^'+^a-/^ 1 )] 


y0 + 2 


y0 + l 


+ ^■r 2 [zi 1 // 2 +^ 2 (i-// 2 )]-K = o 


(10) 


Let T ' be the positive real root of the above equation (10), then T is the optimum cycle time. 

8 2 C(T ) I 


It can also be seen that the sufficient condition for the minimum cost 


> 0 would be satisfied. 


dT 2 1 T=r 

From expression (6), we have the optimum value of Q as follows 

q * = (\-K) ht '[ i + ^r p ] + A 2 r [ i + (id 

p + i p + i 

The optimum average total cost of C(T) is C( T ). 

Some particular cases 
(i). Absence of deterioration : 

If the deterioration of items is switched off i.e. a = 0, then the expressions of C( T* ) and the optimum ordering quantity Q 

become 

c (T*)= ^ + c{\n +^(1-// )]+ ^-r"[/i 1 // 2 +/L 2 (l-// 2 )] 


and Q* = T*[/ \/U +A,(1 — H )] 
where T* = ^ jlK / ll[A^ // 2 + /^(l ~ fP)\ 
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(ii). The deterioration rate is constant i.e. /? = 1. 

In this case, C( T* ) and Q * are given by the following 


c(7")= ( i + ^r l + i 2 r d + |r}] 

which is the cost function of the inventory model obtained by Pal and Mandal[18]. 
and Q * = (A X -A 2 ) juT*[ 1+ °^-T) + A 2 T* [ 1+ ^T*] 

where the real positive value of T* is obtained by the following cubic equation in T 
hot 3 r 3 3.0/1 3\1 (COT + /l) 2 r 5 2 .1/1 2x, „ A 

— t [\/j, + ^ (i — // ) ] + — - — r [x,// +/t,(i -// ) ] - k = 0 


(iii). 77ze demand rate is uniform( A x = /L, j: 

Here the expression for minimum average cost C( T* ) and the optimum order quantity Q are found to be 


C (T*)= — +cl[l+ 

r 4 p+i 


+ y* r 2 + ^Ctp 

2 (P + l)(P + 2) 


T*P j 


and Q* = 1 T*[ l+-^r* A ] 

' /? + l 

where the optimum cycle time T* is the positive real root of the equation 

hap\ T/j+2 + cap\ Tfj+l + r = q 

P + 2 P + 1 2 


(iv). Uniform demand rate and absence of deterioration ( A x = A 2 ; a = 0 ): 

In this case, the expression for the minimum cost and ordering quantity are the following 
C(T ) = c A x + y j2hKA x and Q* = y j2KA x / li 


where the optimum cycle time T is given by 

r= ^2 K/hA, 

Case-2: z is a random point of time. 
From the expression (8), we have 


K c . , a B , a 

C(T, z) = — + - [(^-Aj)z{ 1+ — - Z P } + A,T{ 1+ — — 

r r p+i - p+i 


T p }} 


h 


2 a,p 


+ — [(1-A 0 )z 2 { l + - 
2T (p + l)(p + 2) 


z p } + A 2 T z { 1 + 


2 ap 


(P + l)(P + 2) 

Here the above cost function is a random variable with respect to z. 
Therefore the expected total cost per unit time becomes 


7 ^}] ( 12 ) 
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4<7’) = j + j[{\-X 2 ){E(z) + -^-E(z p+x )} + A 2 T{ 1 + -^T p }] + 


p+\ 


h . . 2 2 aB 

— [(l-A 2 ){E(z )+ — 

2T 1 2 (/? + !)(/? + 2) 


E(z p+2 )} + K T 2 { 1 + 


J3 + 1 
2 a/3 


(/? + !)(/? + 2) 


r A }](i3) 


Let us assume the distribution function of z to be F (■) . 

Therefore we can take 

F (m) = »>j- as in Dhandra and Prasad [17], 

Hence the expected values of z, Z 2 ,Z P+X and z P+ ~ are the following 

E(z) = - £ F Y (u)du ,E( z, 2 ) = - £ 2 uF Y (u)du 
E(z p+ ') =- C (P+l)u p F^(u)du and E(z p+2 ) =- f (P+2)u p+ %(u)du . 

J 0 J 0 

Substituting the above integrals in the equation (13), the expected total cost reduces to the following 

1 p7 — j -t — f r — 

4(T) = — [ K - (A r A 2 ){ c F Y {u)du +h uF Y (u)du +ca w F Y (u)du 


haP |-r 

J+ 


- £ 

\ Jo 


u p+l Fy ( u)du }] + — [ cT + —T 1 + 


A, 


h. 


ca 


T p+X + 


haP 


T P+2 \ (14) 


T ' 2 p + 1 (p + l)(p + 2) 

The equation (14) is of the form 

R + w(T) 

E(T) = , T > 0 where R=K. 

<KT) 

Now (f)(0)) = !//(()) , and (f)(T) , y/(T) are continuously differential functions in [0, oo) with strictly positive derivatives 

(j)\T) and l//'(T) for all T > 0. 

V\T) 

Moreover D(T) = Y 

<f>'(T) 

= [/l 2 + a 2 -A)R(F)][c + hT +ca T p + ^-T p+1 ] 

P + 1 

is a monotonic increasing function of T for A 2 >A l . 

Hence by using Puri and Singh’s[19] result, we get that there exists unique global minimum for a value of T which is the 

optimum cycle time T* . 

Putting T= T in the expression ((14), the minimum expected total cost can be evaluated. 

If A 1 < /lj , we assume the time instant z to be a function of some parameter ‘ V ’ which ranges over some space ‘ F ’ in which 

a probability density function ‘p( y )’ is defined such that J p(y)dy =1 . 

r 

Therefore the expected total cost (f(T) per unit time can be written as follows 


40 T) = [ K+c ( K - ^ 2 ){ J z(y)p(y)dy + J z p+l (y)p(y)dy 

1 r P + 1 i 


} +c/L,T(1h 


a 


P + 1 


T P)1 


2 T 


[(A l -A 2 ){^z 2 (y)p(y)dy- 


2ap 


(P + \)(P + 2) 


\z p+2 (y)p(y)dy 


363 


www.erpublication.org 


Optimal Inventory Management for Items with Weibull Distribution Deterioration under Alternating Demand Rates 


+ /L T 1 { 1 + 


lap 


T P }] 


(P + l)(P + 2) 

From (6), the expected ordering quantity becomes 


05) 


<q> = ( V A, )[ J z(r)p(r)dy + I zP+l (y)p(y) d y ] + a 2 t(i+ -r /; )] ci6> 


If the random variable z (y) and probability density function p( y ) are prescribed, then the right hand side of the equation (15) 

can be evaluated. 

(3 ^ (T^) 

The necessary condition for the minimum expected total cost <^( T ) per unit time is — = 0. 


dT 


This gives the optimum value of T= T* (>0). 

d 2 £(T) 


For 1 = 7 , the sufficient condition for minimum 


dT 2 


> 0 would be satisfied. 


Then the optimum expected ordering quantity of <Q> i.e. < Q* > can be obtained from the expression (16). 


IV. Numerical Examples 

Example- 1 [Deterministic demand pattern i.e. z=known point 
of time] 

Let the parameters of the inventory model be as follows: 

h = $ 3 per unit per year, c = $ 5 per unit , K=$ 20 per order, 
a = 0.002, 

(3 = 2, A l = 10 units, /l 2 = 5 units, p = 0.2. 

Based on these input data, the optimal values of cycle time, 
ordering quantity and minimum total average cost are found 
to be 

T = 1.592 years, Q* = 9.57 units and C( T* ) = $ 55.04 
per year. 

Example-2[Stochastic demand pattern i.e. z=random point of 
time] 

Let z( y ) = a + b V , a, b>0 and the probability density 
function p( y ) be defined as follows 

v(V)= i yexpe/) ' r ~° 

r\/ J 0, elsewhere 

Assuming a = 0.2 and b = 0. 1 and retaining the values of 
the other parameters unchanged(as in example- 1), on solving 

d£(T) 


the equation 


dT 


■ = 0 where ^(T) is given by (15), 


find the optimum mean value of T is 2.035 years. 

From (16) and (15), and using the above mean value of 

T, the optimum expected ordering quantity < Q > and 

minimum expected total cost i ^(T*) are 12.20 units and $ 
55.74 per year. 
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